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ABSTRACT 



Context. The determination of the lower-end of the initial mass function (IMF) provides strong constraints on star formation theories. 
Aims. We report here on a search for isolated planetary-mass objects in the 3 Myr-old star-forming region IC 348. 
Methods. Deep, narrowband CH^o^i and CH^on images were obtained with CFHT/WIRCam over 0.1 1 sq.deg. in the central part of 
IC 348 to identify young T-dwarfs from their 1.6 /jm methane absorption bands. 

Results. We report three faint T-dwarf candidates with C/f4on-C//4off colours >0.4 mag. Extinction was estimated for each candidate 
and lies in the range Aj/ ~ 5 - 12 mag. Comparisons with T-dwarf spectral models, and colour/colour and colour/magnitude diagrams, 
reject two of the three candidates because of their extreme z' - J blueness. The one remaining object is not thought to be a foreground 
field dwarf because of a number density argument and also its strong extinction Ay ~ 12 mag, or thought to be a background field 
T-dwarf which would be expected to be much fainter Models and diagrams give this object a preliminary T6 spectral type. 
Conclusions. With a few Jupiter masses, the young T-dwarf candidate reported here is potentially amongst the youngest, lowest mass 
objects detected in a star-forming region so far. Its frequency is consistent with the extrapolation of current lognormal IMF estimates 
down to the planetary mass domain. 
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To date more than 500 L-dwarfs and about 150 T-dwarf^ have 
been detected using a variety of optical and IR surveys. The vast 
majority of de tected brown dwarfs are old (3-6 Gyr) field dwarfs 
dFahertv et ah, 2009) , having cooled over time. The first detec- 
tion of a T-dwarf in a star-f orming region was put forward by 
IZapatero Oso rio et al.' ( l2002h. named S Orionis 70 and detected 
in cr Orionis. Burgasse r et alj (|2004|) raises questions about the 
authenticity of the object as a cluster member and suggested that 
it was a field dwarf in the line of sigh t, whi ch was backed up with 
a stati stical analysis. iMartfnl (|2004|) and IZapatero Osorio et"an 
( l200 8l) undertook further confirmation of the object by follow- 
up membership and proper motion work. The population of these 
young, low-mass objects is critical to furthering our understand- 
ing of the very low-mass end of the initial mass function (IMF), 
star-forming regions are well suited to search for the lowest 
mass brown dwarfs as young objects are hottest, and so brightest 
immediately after formation. Their temperature decreases over 
time, since by their very nature they are not massive enough to 
begin fusion via hydrogen or deuterium synthesis. 



One ideal system in which to observe young brown dwarfs 
and constrain the IMF is IC 348, a star-forming region to- 
wards the direction of Perseus, centred at (J2000) 03H4™34', 
+32°09'8", and embedded in the foreground part of the Per OB2 
associati on. The age of IC 348 has been determined to be ~ 1-3 
Myr (M uench et al.ll20()3h . This star-forming region is relatively 
nearby and whilst there is controversy surrounding its dis tance, 
IC 348 is located between 261+g pc jScholz etal.lll999 ) and 
316 pc (Herbig 1998) from the Sun whilst the OB Per associa- 
tion , of which IC 348 is a member, i s taken to be between 315 
pc (iLuhman et al.ll2003h and 340 pc ( ICernis|[T993 | )^ An average 
value of 300 (±15 pc) has been used for this work l^rbst 2008]). 
Extinction maps for IC 348, ranging from ~2< Ay<20 mag de - 
pending on the clus ter re gion, were derived by ICernisI (Il993h . 
iMuench et all (l2003h . and lCambresv et al.l (l2006h . 
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iMuench et all ( |2003|) derived the IMF of IC 348 down to 
35 M jiip and found it to be similar to the IMF of the Trapezium 
cluster, having a mode between 0.1-0.2 Mq. Approximately 15- 
25% of the population of the cluster appear to be brown dwarfs 
and their spatial density is independent of the distance from the 
cluster centre. The lowest mass objects were not detected in 
their work because of the detection limits. Those objects with 
a mass of a few M y„p are very faint and cool (~1200 K) enough 
that methane can form in their atmospheres. The presence of 
methane absorption bands in substellar objects defines the L/T 
dwa rf boundary ([B urrows et al. 2006). This was first exploited 
by Main zer & McLeani 0200 3) to look for very faint objects in 
IC 348 to a depth of H ~19.5 mag. Narrow-band 1.65 //m 
methane imaging of IC 348 was conducted in order to filter out 
the lower mass T-dwarfs from the L-dwarfs, by their H - CH4 
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Table 2. Filters and central, effective, wavelengths used in this 
work. 



Filter 






AJAy 




yum 


yum 




z' 


0.88 


0.27 


0.52 


J 


1.25 


0.16 


0.30 


CH^oS 


1.58 


0.10 


0.21 


H 


1.63 


0.29 


0.20 


CH^on 


1.69 


0.10 


0.18 


K, 


2.15 


0.32 


0.13 



colour differences. They found 12-15 candidates which could be 
late M, L or T-dwarfs, with masses > 5 M j„p, although none have 
been confirmed spectroscopically because of their faint nature. 
We report here on a new and deeper methane imaging survey us- 
ing 2 narrow-band filters, C//4on and C//4off, in order to better 
distinguish between T-dwarfs and more massive young stellar 
objects in IC 348. 

In Section 2 there appears a short description of the obser- 
vations and instrumentation used, and of the methane filters. We 
also describe the photometric extraction method and the result- 
ing catalogues. The criteria for the selection of the candidates are 
presented in Section 3, and the promising candidates identified. 
Section 4 discusses the photometric properties of the candidates, 
their nature and membership status, and Section 5 follows with 
conclusions. 



2. Observations and photometric extraction 

This work has made use of the 3.6m Canada France Hawaii 
Telescop^ (CFHT) on Mauna Kea, specifically the infrared 
camera WIRCam (PugeLet al. 2004) and the wide field imager 
MegaCam (Boulade et all i2003i) . Basic information regarding 
these surveys is summarised in Table [T] whilst the central wave- 
lengths for the filters used as well as the corresponding extinc- 
tion coefficients can be found in Table|2] Central, effective, wave- 
lengths and extinction coefficients are taken from the Spanish 
Virtual ObservatorjQ. 

2.1. WIRCam 

WIRCam observations were obtained under programme 06BF23 
in service mode between September 2006 and January 2007. 
A grid of four IR arrays of 2048x2048 pixels make up the 
WIRCam detector yielding a field of view of 20 'x20 '. The 
pixel scale is 0.306 " or 0.15 " with microdithering (used for 
the J, H and Kg images). WIRCam was used to take narrowband 
C//4off and C//4on along with J, H, and Ks images of IC 348. 
The C//4on and C//4off pointings are both centred at the same 
location, which includes the cluster's center, while four over- 
lapping WIRCam fields were taken in J, H and to provide 
a larger areal coverage (see Table [1]). Total integration times of 
3.7h in C//4ofF and 1.4h in C//4on were obtained from multiple 
individual exposures of 30s, acquired using a 7-position dither- 
ing pattern to fill in the 45"-wide gaps in between IR arrays. The 
same observing strategy was used for JHK, images, with indi- 
vidual integration times of 45, 10, and 15s, respectively, yielding 
total integration times of 1260s (J), 560s (H) and 420s (K^) for 
each pointing. 

2 www.cfht.hawaii.edu 

^ http://svo.laeff.inta.es/theory/filters/index.php 



Individual images were detrended using the 'i'iwi pipeline 
(Albert et al., in prep.) at CFHT and sky subtraction, stack- 
ing, photometric and astrometri c calibrations , and quality con- 
trol were performed at Terapi?iPI (lMarmol2007b . The seeing (PSF 
FWHM) measured on the final images was between 0.55" and 
0.65". The co-added JHKj images were photometrically cali- 
brated using the 2MASS catalogue over the same area and renor- 
malized to an arbitrary photometric zero-point of 30 mag. The 
methane images have no external photometric calibration and 
the CH4 magnitudes are given here on an arbitrary albeit inter- 
nally consistent scale, so that C//4on-C//4off^0 for unreddened 
field dwarfs (see below). The Terapix pipeline additionally pro- 
duced a catalogue of objects with aperture photometry for each 
filter 

2.2. MegaCam 

MegaCam observations were obtained under programme 
06BF28 in queue service mode on September 21-23, 2006. 
MegaCam consists of a grid of 36 2kx4k CCDs covering a 1 ° x 1 ° 
footprint. The pixel size is 13.5 fim and the pixel scale is 0.185". 
MegaCam was used to take 6 dithered z'-band images of 1500s 
each, yielding a total integration time of 9000 s. After detrending 
the images at CFHT, sky-subtraction, stacking, and photometric 
and astrometric calibrations were done at Terapix. The images 
were photometrically calibrated using standard stars routinely 
observed by the Queue Service Observing team at CFHT. The 
seeing during the observations ranged from 0.65" to 0.80". 

2.3. Photometric catalogues 

A combination of SExtractor dBertin & ArnoutsI Il996h and 
PSFex (Bertin et al., in prep.) was used on each image to ex- 
tract sources and build a photometric catalogue for each image. 
On a first step, SExtractor extracts well defined stellar-like ob- 
jects, which are used by PSFex to compute a PSF model that is 
allowed to vary with position on the detector. Then, SExtractor 
uses this PSF model to more accurately extract and measure the 
photometry of all the sources detected on the image. Objects 
were detected at 3.5 cr above the background noise (see Section 
2.4 for further details), most notably to avoid the high incidence 
of nebulous emission and reflected nebulosity in the region. This 
threshold is reasonably low because our interest lies in detect- 
ing the faintest possible objects in the region, the methane T- 
dwarfs. Unfortunately, this low threshold also increases the inci- 
dence of non-stellar-like detections because of nebulosities, and 
many misdetections created via crosstalk in the detector (though 
this was only particularly problematic in the narrowband obser- 
vations). 

Further analysis was conducted using publicly-available 
2MASS data for those detected objects that were bright enough 
to have a 2MASS counterpart. This was useful for two rea- 
sons: firstly that the photometric accuracy and pipeline reduc- 
tion of the WIRCam images could be compared with those from 
2MASS; and secondly that the difference in the photometric sys- 
tems could be estimated from the photometric differences be- 
tween the CFHT and 2MASS. The extracted J, H and K, cata- 
logues were matched with the J, H and 2MASS catalogues 
resulting in 206 stars that had overlapping magnitudes. Based 
on these, the dispersion between the extracted J, H and Ks cat- 
alogues and the 2MASS /, H and Ks catalogues was calculated. 
Good agreement was shown by mean magnitude differences of 

terapix.iap.fr 
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Table 1. Instruments used with field of view (FOV), filters, and the central coordinates of the pointings. The 7, H and Ks final 
WIRCam images are made up of four overlapping fields (A, B, C and D). 



1 CiCoHJUC 


FOV 






Trnfir»n tiTni^ (hi 
'IdLlUll LiiiiC yii) 


Pr^intiTifT r'r\r^t"ri 1 n Q tf* t; ( 1/(111(11 
± UlllLillg L-UUl U.llId.lCS yj Z.\J\J\J ) 


CFHT MegaCam 


0.96°x0.94° 


z' 


2.5 




03'■44"'36^00 +32°0l'5a'0 


CFHT WIRCam 


20'x20' 




3.7, 


1.4 


03M4"' 14^80 +32°05'06'.'0 


CFHT WIRCam 


20 'x20 ' 


J, H, K, 


0.35 


, 0.16, 0.12 


03''44"'14^80 +32°05'06'.'0 




A 










B 








03M4"'14^80 +32°15'06'.'0 




C 








03M5"'02^00 +32°05'06'.'0 




D 








03M5™02^00 +32° 15' 06'.'0 



0.08 + 0.02, 0.03 + 0.05, and 0.04 + 0.02 for J, H and K, re- 
spectively. These values hold for the brightest stars in our im- 
ages, since fainter ones are undetected in the 2MASS survey. In 
turn, this means that colour effects arising from CFHT filters are 
not taken into account in this calibration. Indeed, while the pho- 
tometric zero-point used here is obtained from 2MASS, colour 
effects are not corrected for and our WIRCam J, H, Ks, CH non 
and C//4off photometry is given in the CFHT Vega system and 
MegaCam z' in the AB system. 

Each waveband had a catalogue created which was then col- 
lated and cleaned of saturated objects, obvious artefacts, and 
a fraction of extended sources identified from their large PSF 
FWHM. An approximation of the completeness of the photo- 
metric catalogues was estimated from log(Noiy) vs. magnitude 
histograms, where No^y is the number of stellar-like objects de- 
tected on the images. We thus derived completeness limits from 
Figure [2 f top j of -23.5 (z'), 21.5 (7), 20.0 (//), 18.9 (^T,), 20.3 
(C//4on), and 20.7 (C//4ofF). SExtractor photometric errors for 
each of the six bands are plotted in Figure[T]ftoftomj, highlight- 
ing the increasing photometric error for the faintening of the ob- 
jects' brightness. Uniquely, the CHtioS error terminates at 0.06 
magnitudes as the C//4off image was used for the detection of 
objects over all bands (see Section 2.4 for further details). 



2.4. Methane photometry 

The presence of methane absorption bands in the near-IR spec- 
trum of T-dwarfs can be used to identi fy T-dwarf candidates pho- 
tometrically (e.g. iTinnev et al.ll2005l) . Narrowband and 
CHuon data was taken in order to classify these objects. CH40S 
measures the pseudo-continuum at 1.58 /im while CHuon sam- 
ples the methane absorption band at 1.69 fiva. The passband of 
the two filters is overlain onto the spectra of a TO. 5 and a T8 
dwarfs in Figure |2] Note the greater methane absorption in the 
T8 dwarf spectrum compared to the TO. 5 spectrum in the region 
around 1 .69 //m. 

Spectrj^ of field dwarfs from LI to T8 were convolved 
jDelorme et al.l l2008') with the WIRCam CH^on and CHaoS fil- 
ters and the resulting methane colours plotted against spectral 
type in Figure |2] The methane colours are seen to smoothly 
increase towards later spectral types. Whilst L-dwarfs have 
C/f4on-Ci/4off colours equal to zero, T-dwarfs have colours 
above 0.1 mag which rapidly increase towards later T-types. 
To date, no other types of objects in current knowledge has 
this sequence of methane colours. Thus, the C//4on-C//4off 
colour provides a useful means to separate L and T-dwarfs. 



Additionally, as late-type T-dwarfs may remain undetected in 
the CH/^ovi image, and due to strong CH4 absorption, we first 
performed object detection at 3.5cr on the CHiioS image, then 
performed PSF photometry at the location of these objects on 
both the and CHuon images. This ensures that all ob- 

jects that could possibly be T-dwarfs are searched for in both 
filters. 

Figure |2] provides an empirical calibration of the methane 
colours against spectral type for field dwarfs. For comparison. 
Figure |3] shows the CH^on-CHnoS colour as a function of ef- 
fective temperature {Te fr) as predicted by COND and DUSTY 3 
Myr and 5 Gyr models (Ichabrier et al.ll2000HAllard et aLllMTt 
Baraffe et al. 2003). While the DUSTY models are roughly sim- 
ilar at both ages, with C//4on-C//4off ~ at T^fj > 1500 K, 
the COND models predict bluer colours for younger objects at 
Tef / < 1 500 K. According to the models, IC 348 T-dwarfs would 
then have a smaller T^jf and thus a later spectral type than field 
T-dwarfs for the same C//4on-C//4off colour. 

3. Results 

Catalogues were obtained as described above for each of the 
MegaCam/WIRCam images. The COND and DUSTY 3 Myr 
and 5 Gyr model^ are also in the CFHT Vega system where 
the z' band has been converted into AB magnitudes for our pur- 
poses. Similarly, the empirical field L-T dwarf sequence pho- 
tometry has been adjusted to be consistant with the models and 
our data. In the colour/magnitude diagram presented the field 
dwarf sequence has been shifted to the distance of IC 348, taken 
to be at 300 pc. 

3.1. T-dwarf candidate selection 

Figure |4] shows C//4on-C//4off against C//4off for stellar-like 
objects detected on WIRCam methane images. An object was 
considered as a T-dwarf candidate if C//4on-C//4ofF > 0.4 
mag, which corresponds to > 3.5 cr above the L/T tran- 
sition (C//4on-C//4ofF - 0) for the faintest detected ob- 
jects. Assuming that young T-dwarfs follow the field dwarf 
methane relation, the spectral type of those objects with 
C//4on-C//4ofF = 0.4 mag is close to T3 from Figure |2] 

Objects with C//4on-C//4off > 0.4 mag were initially all 
classified as T-dwarf candidates. The objects just below this limit 
could still be T-dwarfs from their C//4on-C//4off colours as 
shown in Figure |2] but would be more difficult to extract be- 
cause of photometric errors. With this criterion 136 sources were 
selected. This sample was then visually scrutinised on both the 



5 www.jach.hawaii.edu/~skl/LTdata.html 



http://phoenix.ens-lyon.fr/simulator/index.faces 



4 



A. S. M. Burgess et al.: Young T-Dwarf Candidates in IC 348 




2 - 



18 20 22 24 26 16 18 20 
z / mag CH^off / mog 



16 18 20 22 
CH^on / mag 



0.1 
0.08 

oi 

£ 0.06 
\ 

o 0-04 
oj 

0.02 


0.1 
0.08 

en 

£ 0.06 
o 0.04 
0.02 



















- ■■■> i - 
: f- '- 
















- :W - 


^, , , l' , , , 1 , , , 1 , 





" 


18 20 22 


18 20 22 


16 18 20 


J / mag 


H / mag 


Kg / mag 









" T. . \ . . . \ . . . 







20 22 
z' / mag 



16 18 20 
CH^off / mag 



16 18 20 22 
CH^on / mag 



Fig. 1. top: log(Nobj) vs magnitude (where Nobj is the num- 
ber of detected objects in each band): The completeness lim- 
its are found from these plots for the six imaged bands, z, J, 
H, K, CH4on and C//4ofF. The line has been fitted to the his- 
togram to find where the turning point occurs, thereby indicating 
our completeness limit, bottom: Photometric error, as measured 
by SExtractor, against magnitude. Note that the J, H, Ks images 
have multiple tracks as they are composed from 4 fields which 
were not acquired at the same time. 



CH40S and C//4on images before being reduced to a shortlist 
of 12 possible candidates. Upon visual inspection, we found that 
the rejected 124 'objects' were saturated stars, nebulous detec- 
tions or ghosts. Of these 12 possible candidates, the shape of the 
PSF and the contours were further examined using IRAF, where 
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Fig. 2. C//4on-C//4ofF colour vs. spectral type. The triangles 
are computed methane colours for a spectral sequence of field 
L and T dwarfs (see text). The dotted line at C//4on-C//4ofF= 
0.1 mag highlights the limit where methane absorption becomes 
conspicuous in the photometric index, which occurs at a spectral 
type TO. Note how rapidly the methane index increases towards 
later spectral types. Inset: Dwarf TO. 5 and T8 spectra are shown 
in the region of the 1 .6 /jm methane absorption band. The pass- 
band of the WIRCam C//4off and CH4on filters is superimposed 
(dashed). 



9 were identified as ghosts or detector cross-talk, so reducing 
this figure to 3 likely candidates. Thumbsized CH40S images of 
the three candidates are shown in Figure |5] 



The PSF photometry of the 3 candidates in z', J, H, 
Ks and C//4on/C//4off filters is listed in Table |3] Candidate 
IC348_CH4_2 was not detected in the z'-band. The z'-band de- 
tection limit was estimated by simulating a set of 19 48 stars wit h 
magnitudes between 20 and 26 using Skymaker (iBertinll2008l) . 
These synthetic stars were then stacked o nto the CH40 S foot- 
print of the original z' image using Swarp (Marmo 2007). At the 
3cr level, 1688 objects were detected where the faintest had an 
input magnitude of 25.7, which we take as being the detection 
limit in the z' filter. 

Most of the objects detected on the WIRCam images lie 
along the line C//4on-Cj/4off - in Figure |4] as expected 
for field dwarfs jTinnev et al.ll2005l) . the photometric error in- 
creasing with magnitude. The 3 candidates are located in the 
faintest region of the plot, whilst also having increasing methane 
colours. The C//4on-C//4ofF rms photometric error at the mag- 
nitude of the candidates is cr ~ 0.12 mag. The 3 candidates have 
C//4on-C//4off colours of 0.54, 0.78, and 0.49 mag, respec- 
tively, corresponding to a detection level of > 4cr (see TablelDl. 

Note that a number of objects in Figure|4]have strongly neg- 
ative CH4on-CH4off colours, well beyond the rms photometric 
error Upon inspection of the images, all these objects turn out to 
be young IC 348 members deeply embedded in bright compact 
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Fig. 5. Thumbsized Ci/4off images of the three T-dwarf candidates, left to right, IC348_CH4_1, IC348_CH4_2 and IC348_CH4_3, 
respectively. Black circles highlight the location of the T-dwarf candidates. The images are 48" to a side, where North is up and 
East is to the left. 

Table 3. PSF photometry and photometric errors of the three T-dwarf candidates, in magnitudes. 
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z' 
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o-j 


H 
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0.07 


21.62 


0.04 
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0.08 
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Fig. 3. C//4on-C//4ofr colour vs. Tejf. 3 Myr COND {thick 
solid), 5 Gyr COND (thin solid) 3 Myr DUSTY (thick dashed) 
and 5 Gyr DUSTY ( thin dashed) models are shown. The COND 
models predict later spectral types for 3 Myr IC 348 objects than 
for field dwarfs for given methane colours. 



nebulosities. The large extinction they suffer results in strong 
reddening, thus yielding blue (i.e. negative) C//4on-C//4off 
colours in spite of the small wavelength difference between the 
2 narrow-band methane filters. Reddening also accounts for the 
asymmetric colour distribution of the faint background objects 
in this plot. 

Finally, a cross-compa rison was made with the study of 
iMainzer & McLeanI (l2003l) . They classified 5 candidates as be- 



en 
o 



0.5 



O 
I 

c 
o 

X 

o 




-0.5 



CH^off / mag 

Fig. 4. CH4on-CH40« colour vs C//4ofr. The large (blue) 
squares show the three T-dwarf candidates with C//4on-C//4off 
> 0.4. The extinction vector is shown for Ay = 10 mag. The 
sources with strongly negative colours are deeply embedded ob- 
jects (see text). 



ing possibly M, L or T-dwarfs, all of which are detected in 
our images but with C//4on-C//4off colours between -0.14 and 
0.04 mag, placing them out of our T-dwarf candidate criteria. 



6 



A. S. M. Burgess et al.: Young T-Dwarf Candidates in IC 348 



3.2. Reddening and spectral type estimates 

Estimation of the candidates' extinction is required in order to 
be able to calculate their absolute magnitude and estimate their 
spectral type. Extinction has been estimated using colour/colour 
diagrams of C//4on-C//4off versus J - H, J - Ks and H - Ks, as 
plotted in Figure|6l The extinction was computed for each candi- 
date using the extinction vector and regressing the objects back 
towards the 3 Myr COND model. The final extinction value is 
the average of the 3 results obtained from each colour/colour di- 
agram and is summarised in TableH] along with C//4on-C//4off 
colour, detection level and spectral type for the three candidates. 
The estimated extinction is an upper limit if the candidates be- 
long to the cluster as their true dereddened colours are probably 
intermediate between the 3 Myr COND and DUSTY models. If 
they are field dwarfs however, the given value is a lower limit as 
the objects should be dereddened towards the 5 Gyr field dwarf 
sequence that is bluer than the 3 Myr COND model. 

We then used these results to compute the dereddened 
C//4on-C//4off colour for each candidate. We found these val- 
ues to be 0.69 + 0.16, 1.15 ± 0.24 and 0.76 + 0.15 mag 
for IC348.CH4.1, IC348_CH4.2 and IC348.CH4_3, respec- 
tively. From Figure |2] these dereddened colours correspond to 
a spectral type of T5+° ^ T6+i 5, and T5+[5 ;^ for IC348_CH4.1, 
IC348_CH4_2 and IC348.CH4_3, respectively. As noted above, 
this may be a lower estimate as the models suggest that young 
T-dwarfs have a lower effective temperature than field T-dwarfs 
for the same C//4on-C//4off colour (see Figl3]i. Moreover, if the 
candidates were 5 Gyr field objects, they would have a larger ex- 
tinction and therefore an even larger intrinsic C//4on-C//4off 
colour, which would also yield a later spectral type. 

4. Discussion 

We discuss here the likelihood that the methane candidates re- 
ported above are bona fide young, very low-mass members of 
the IC 348 star-forming region instead of being more evolved 
field T-dwarfs located on the line of sight to the young cluster. 
We also compare the number of T-dwarf candidates we identify 
in our survey to the number of expected planetary mass objects 
in IC 348, by extrapolating recent estimates of the substellar IMF 
to the planetary mass domain. 

4.1. Membership 

Further colour/colour and colour/magnitude diagrams (CMD) 
were plotted to constrain the candidates' status. Figure [7] shows 
the J/J-Ks CMD. Also plotted are the synthetic magnitudes and 
colours of the field T-dwarfs when shifted to the cluster distance. 
The 5 Gyr DUSTY and COND models are drawn to highlight 
the agreement between the models and the field dwarf sequence. 
The field T-dwarf sequence begins in between the 5 Gyr DUSTY 
and COND models, before sweeping towards the 5 Gyr COND 
model at ~1500 K that sharply increases in faintness with the 
late-T set. 

In this diagram the three candidates appear to be confused 
with the late-L, early-T field dwarfs. However, once they are 
dereddened using the extinction values given in Table [H they 
become bluer and brighter All three candidates appear to follow 
the 3 Myr COND model quite closely. The fact that the dered- 
dened candidates appear brighter than field T-dwarfs shifted to 
the cluster distance suggest that they are younger indeed. A 
young T-dwarf has a larger radius than a field T-dwarf for a given 
spectral type as it is still contracting. According to the COND 



■ 

■ 




J-K^ / mag 

Fig. 7. J vs J - Ks CMD. The (blue) squai-es are the three T- 
dwarf candidates as observed, whilst the (magenta) pentagons 
represent their dereddened location in the CMD. Field T-dwarfs 
(red triangles) as well as S Ori 70 (black square, labelled) shifted 
to the distance of IC 348 are shown for compariso n. Dereddened 
M6-M9 IC 348 dwarfs froml Luhman et al.1 (12003 ^ are plotted for 
comparison (small green squares). DUSTY and COND models 
are shown as in Fig. [3] 



models, the difference amounts to about a factor of 2 in radius 
between 3 Myr and 5 Gyr (for T,//= 1000- 1500 K), which re- 
sults in a 1.5 magnitude increase, consistent with the observed 
location of the candidates in Figure|2l 

There is a possibility that the candidates are in fact field 
T-dwarfs located at a closer distance (between ~100 and 200 
pc) along the line of sight to IC 348. The probability of one 
of the candidates being a field T-dwarf instead of an IC 348 
cluster member can be estimated independently from extinction. 
According to Figure |2] the spectral type range T3-T5.5 corre- 
sponds to the measured C//4on-C//4off colours of our candi- 
dates at an age of 5 Gyr. Thusly, using the number density of 
T3-T5.5 field dwa rfs in the solar neighbourhood, i.e. ~1 per 714 
pc^ according to Metchev et al.l (2008) and the footprint of the 
CH4 image of 0.1 1 sq.deg. gives an estimate of 0.1 1 ± 0.06 T3- 
T5.5 foreground field dwarfs in the direction of IC 348. As we 
are concerned here with the probability of one of our candidates 
being a foreground field dwarf, this is a fairly robust method for 
determining the population density for field dwarfs in the direc- 
tion of IC 348. However, this estimate is put into context when 
taking into account the large extinction values of the three can- 
didates estimated from Figure |6l Even the least extincted candi- 
date, IC348_CH4_1, has four magnitudes more extinction than 
expected for a foreground field dwarf (<1 mag), so all three ob- 
jects must be near to or behind IC 348. However, the candidates 
cannot be background field T-dwarfs, seen through the IC 348 
cloud, as their luminosity would then be much too high for their 
estimated spectral type. Finally, as indicated in Table|4] all three 
candidates are located within the cluster's boundary (4' core ra- 
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Fig, 6. C//4on-C//4off \s. J - H, J - Ks and H - Kg. T-dwarf candidates are plotted as (blue) squares. Photometric error bars are 
from SExtractor. Field T-dwarfs (red triangles) are shown for comparison. Lines show the 5 Gyr COND and DUSTY, and 3 Myr 
COND and DUSTY models as in Fig.|3] Candidates were dereddened towards the 3 Myr COND model using the extinction vector. 



Table 4. Summai-y of values for the three candidates, where the cluster centre is taken to be at 03''44'"34" +32°09'48'.'0 (J2000). 



lAU name 


C//4(on-ofF) 


Det. 


Ai//mag 


Est. sp. 


Distance from 


Object coordinates 




(mag) 


level 




type 


cluster centre 




CFHTJ0344+3202_(IC348.CH4.1) 


0.54 


4.3 cr 


5.0 ± 1.2 


T5!;'-' 


7'35" 


03H4'"49^24 +32°02'48'.'4 


CFHTJ0344+3206_(IC348.CH4.2) 


0.78 


6.2 or 


12.4 ± 3.9 


T6+' 

^ "-0.5 


4' 20" 


03''44"'49^52 +32°06'35'.'4 


CFHTJ0344+3156_(IC348.CH4.3) 


0.49 


A3 a 


9.0 ± 1.3 


T5 ± -0.5 


13' 56" 


03h44"57\95 +31=56' 43'.'3 



dius, 10-15' halo; Herbig 1998, Herbst 2008) and are thus spa- 
tially consistent with being IC 348 members. 

4.2. Contaminants 

A further, useful diagram in defining these candidate objects is 
the z' - J/ J - H colour/colour diagram shown in Figure[8] (note 
that the z' is in the AB system). Here it can be seen that the two 
detected dereddened candidates are much bluer in z'-J than both 
the field dwarf sequence and the COND 3Myr and 5Gyr mod- 
els. There is mounting empirical evidence for young T-dwarfs to 
have bluer z'-J colours than field T-dwarfs because of the effects 
reduced gravity has on the opacities (P. Delorme, priv. comm.). 
This effect stems from the strong potassium KI (7687&7701 /vm) 
doublets, whose wings fall within the z' band. The lower the 
gravity the lesser the line broadening of these elements so less 
flux is lost by absorption in the z' band, resulting in the ob- 
served bluer z'-J colours (F. AUard, priv. comm). However, 
these effects are unlikely to explain the extreme blueness of two 
of the three objects and suggests that candidates IC348_CH4_1 
and IC348_CH4_3 are unlikely to be T-dwarfs. The final remain- 
ing candidate, IC348_CH4_2, remains a good T-dwaif candidate 
because of its non-detection in z' ■ 

If candidates IC348.CH4_1 and IC348.CH4_3 are not young 
T-dwarfs then what are they? We checked whether extragalactic 
objects could contaminate this region of the diagrams. To this 
end, galaxie^l from SWIRE, 2MASS and SDSS were found to 
have very different colours (z' - / ~ 6-8mag) to our candidates. 
Similarly, quasars (Labitaet al. 2008) also have J-band magni- 
tudes of ~ 24-25 mag, making them much fainter than our ob- 
jects. Late-type (emission line) galaxies appear tightly in the re- 
gion J-K ~ 0.95 mag and J-H ~ 0.88 mag and so again cannot 
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Fig. 8.z' -JvsJ-H diagram. DUSTY and COND models are 
shown as in Fig. [3] Symbols shown as per Fig.|2] 

be confused with our candidates ("Chang et al."2006'). The status 
of the two rejected candidates continues to be unclear. 

4.3. Comparison of IC348.CH4^ with S Ori 70 

IZapatero Osorio et al.l ( l2002t) found S Ori 70 to be a faint mid-T 
type object towards the direction of the young cr Orionis clus- 
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ter, with an estimated mass in the interval 2-7 My„p. However, 
there has been some contention with regard to the member- 
ship of S Ori 70 to th e cr O rionis cluster. Further work by 
IScholz & Javawardhanal (1200 8l) indicated that this was a true 
cluster member by looking at the IRAC data, and suggested that 
there is a possibility of a disk structure which gives rise to the 
MIR excess. 

The J/J - K, CMD in Figure [T] has the IC 348 dereddened 
candidates plotted alongside S Ori 70 for c omparison. The J, H 
and K s data for S Ori 70 was taken from Zapatero Osorio et alj 
( |2008|) and adjusted to the CFHT photometric system using the 
3 Myr COND CFHT and 2MASS models. The differences for 
'^iMASs - ^CFHT due to colour effects are 0.6 mag, -0.1 mag, 
and -0.05 mag, for J, H and Ks respectively. In this diagram, 
candidate IC348_CH4_2, for which we estimated a spectral type 
of T6^y J, falls at approximately the same location as S Ori 70, 
shifted to IC 348's distance. 

Figure|9]shows the J - H/H - Ks colour/colour diagram. The 
dereddened J — H colours of remaining IC348_CH4_2 T-dwarf 
candidat e is -0.44 which cor responds to a spectral type of around 
T6-T8 (Luh man et al.ll20d7l) and so is consistent with our previ- 
ous estimate. Again in this plot, the location of IC348_CH4_2 is 
similar to that of S Ori 70. A remarkable property of these two 
objects is that they have H - Ks colours too red for their esti- 
mated spectral types, yet also have blue J — H colours consistent 
with them being T6 dwarfs. According to^ Zapatero Osorio et alj 
(i2008.) . this may be a signature of youth. 

4.3.1. Spitzerdata 

We searched the Spitzer public archive for IRAC data. IC 348 
was the prime target of programmes 6 (c2d, RI. Fazio) and 
36 (RI. Fazio). The c2d data are made up by shorter expo- 
sures, resulting in our use of the deeper (16x1 00s) images of 
programme 36. We retrieved the data and processed them us- 
ing standard procedures with the recommended MOPEX soft- 
ware. All three candidates are detected in the [3.6] and [4.5] 



bands, except IC348_CH4_1 which falls out of the [4.5] band 
field of view. PSF photometry was performed using Starfinder 
(Diolaiti et al. 2000) and the fluxes were translated into magni- 
tudes using the zeropoint fluxes provided by the Spitzer Science 
CenteiQ Measurement uncertainties were tentatively estimated 
from the Poisson noise weighted by the coverage maps of the 
mosaics. The final photometry is given in Table |5] and the final 
errors include both measurement and zeropoint flux uncertain- 
ties. 

In Figure [To] the Ks - [3.6] vs [3.6] - [4.5] colour/colour 
diagram is plotted with the two can didates, S Ori 70 a nd the 
IC 348 M6-M9 dwarf sequence from Luhman et al. (2005). The 
M, L and T IRAC field dwarf sequence from Patten et al. (2006|) 
have also been plotted for comparison. The [3.6] - [4.5] colour 
could be a good indicator of effective temperature, or spectral 
type, giving IC348_CH4_2 a spectral type of ~T5, in line with 
our other estimates of T6. The Ks - [3.6] colour, however, relates 
to a spectral type of LI, which is far too early for such a low- 
mass and cool object. 

TheKs - [3.6] colour of IC348_CH4_2 appears significantly 
bluer compared to the field sequence , which may be because of 
reduced gravity dLeggett et al.ll2007l) . For mid-type and later T- 
dwarfs the reduced pressure broadening of H2, as a gravity ef- 
fect, makes the Ks band brighter, whilst the [3.6] band faintens 
because of additional CH4 absorption at ~ 3 jjm, resulting in a 
bluer Ks - [3.6] colour Another consequence of lower gravity is 
the faintening eff'ect the CO abundances have on the [4.5] band 
for early T-types, thus balancing the [3.6] - [4.5] colour. These 
colour effects of low gravity are in agreement with this object 
being a young IC 348 T-dwarf. 

Other effects such as strong sedimentation and httle or no 
vertical mixing in the atmospher e could also contrib ute to the 
observed bluer Ks - [3.6] colours dLeggett et al.ll2007l) . 

The possible existence of a hot inner dusty disc around S Ori 
70, in keeping with other young members of S Orionis, can give 
rise to its redder [3.6] - [4.5] colour. 

4.4. The lower end of the IMF 



iLuhman et all (l2003l) derived a nearly complete IMF for IC 348 
down to 0.03 M0 for Ay < 4 mag. They find a ratio of brown 
dwarfs (BDs, 0.02-0.08 Mq) to stars of about 12%. lMuench et all 
(I2003h similarily derived an IMF for IC 348 down to 0.04 Mq, 
with a mode between 0.08 and 0.2 M©, finding a ratio of BDs 
to stars of about 14% . From a deep J-band survey of the cluster, 
iPreibisch etlD (l2003h derived a BD to star ratio of 10%. 

According to the COND and DUSTY models, the mass 
range of 3 Myr old T-dwarfs is between 0.001 M0 and 
0.005 Mq, corresponding to masses from ~1 My„p to ~5 My„p. 
However, the models are somewhat uncertain at very low 
masses and young ages. Currently the models can have up 
to a 50% discrepancy in the masses for objects younger than 
5 Myr (Chabrier, priv. comm). Still, the T-dwarf candidate 
IC348_CH4_2 is likely to be less massive than 10 My„p if be- 
longing to IC 348. 

In order to obtain an estimate of the number of objects within 
this mass range that are expected to be in IC 348, we extrapo- 
lated current estimates of the IMF to the planetary mass regime. 
A lognormal esti mate of the field IMF for unresolved systems 
was provided by IChabrierl (|2003|) . where the mode of the dis- 
tribution, mo, is 0.22 Mq and its width, cr, is 0.57. This log- 
normal IMF would predict ~1% objects in the mass range be- 
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Table 5. Spitzer PSF photometry and photometric errors of the three T-dwarf candidates. IC348_CH4_1 appeared outside the FOV 
in the [4.5];um image. 



Object 


F[3.6]/mJy 


F[3.6]_error 


F[4.5]/mJy 


F[4.5]_error 


[3.6]/mag 


[3.6]_error 


[4.5]/mag 


[4.5]_error 


IC348_CH4.1 


0.0035 


0.0001 






19.75 


0.03 






IC348_CH4.2 


0.0093 


0.0003 


0.0136 


0.0004 


18.70 


0.03 


17.80 


0.03 


IC348_CH4.3 


0.0245 


0.0003 


0.0365 


0.0007 


17.65 


0.02 


16.73 


0.02 
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Fig. 10. Ks - [3 .61 y s [3.6] -[4.5]. The field dwarf sequence from 
iPatten et"aLl ( 20061) has been plotted (red crosses) for compari- 
son. IC 348 M6-M9 dwarfs with disks (filled gree n squares) and 
those without (hollow green squares) are plotted (Luhm an et alj 
l2005h . IC348.CH4_2 and S Ori 70 appear significantly bluer in 
the Ks - [3.6] colour than for field dwarfs. The extinction vec- 
tor has also been plotted using values from the Spanish Virtual 
Observatory^. 



low 10 Mjup. Similarily. iMoraux et alJ (l2007h suggested a uni- 
versal lognormal IMF for systems in young open clusters, with 
wo=0.30 + 0.05 M0 and cr=0.55 + 0.03, which would predict 
~0.4% of objects in the 1-10 My„p mass range. iMuench et al] 
j2007^ conducted a Spitzer census of IC 348 and stated that the 
population of IC 348 is in excess of 400 members when taking 
into account unseen diskless members. Based on the above IMF 
and population estimates, we would thus expect about 1 .6-4 ob- 
jects in the mass range 1-10 My„p in the cluster The discovery 
of one T-dwarf candidate close to the completeness limit is thus 
consistent with the extrapolation of current lognormal IMF esti- 
mates down to the planetary mass domain. 

Another estimate can be found using the power law derived 
bylCaballero et al. (20^ and|Lodieu et al. (2009) for cr Orionis 
and lLodieu et all (120(51 for Upper Sco, ^ cc M where a 
= 0.6 for the mass range 0.3 - 0.01 Mq. Using the IMF from 
the Ay < 4 mag selected census of IC 348 by iLuhman et alj 



( l2003 f) and extrapolating the power law for the last two data 
points above their completeness limit towards the lowest mass 
domain gives a different estimate. This corresponds to a much 
larger number of predicted objects than for the lognormal esti- 



mate, 25 + 16 objects in the 1-10 My„p mass range and 15 + 9 
objects in the 1-5 My„p mass range. Even though our survey is 
complete down to Ay ~ 12 mag for T3 - T5.5 dwarfs, we did not 
detect any such object with an Ay < 4 mag. This suggests that 
the current power law approximation over-estimates the number 
of low-mass objects in IC 348. 



5. Conclusions 

From a deep methane imaging survey of the star-forming re- 
gion IC 348 we identified 3 T-dwarf candidates over the area 
of the cluster. After colour/colour and colour/magnitude dia- 
gram analysis two candidates have been rejected for being too 
bright at optical wavelengths. The remaining candidate, has an 
estimated spectral type of T6 and theoretical models suggest a 
mass of a few M /„y, for this object at 3 Myr From its luminosity, 
colour, extinction and spatial location, IC348_CH4_2 is a prob- 
able IC 348 T-dwarf member, and so is among the lowest mass 
objects observed so far in a star-forming region. The frequency 
of isolated planetary mass objects reported here for IC 348 is 
consistent with the extrapolation of current lognormal IMF esti- 
mates to the planetary mass domain. 
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